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Abstract 

We analize the use of algorithms based in x-space for the solution of renormalization 
group equations of DGLAP-type and test their consistency by studying bounds among par- 
tons distributions - in our specific case Soffer's inequality and the perturbative behaviour of 
the nucleon tensor charge - to next-to-leading order in QCD . A discussion of the perturba- 
tive resummation implicit in these expansions using Mellin moments is included. We also 
comment on the (kinetic) proof of positivity of the evolution of hi, using a kinetic analogy 
and illustrate the extension of the algorithm to the evolution of generalized parton distri- 
butions. We prove positivity of the non-forward evolution in a special case and illustrate a 
Fokker-Planck approximation to it. 



One of the most fascinating aspects of the structure of the nucleon is the study of the distribution 
of spin among its constituents, a topic of remarkable conceptual complexity which has gained 
a lot of attention in recent years. This study is entirely based on the classification and on the 
phenomenological modeling of all the leading-twist parton distributions, used as building blocks 
for further investigations in hadronic physics. 

There are various theoretical ways to gather information on these non-local matrix elements. 
One among the various possibilities is to discover sum rules connecting moments of these distribu- 
tions to other fundamental observables. Another possibility is to discover bounds - or inequalities 
- among them and use these results in the process of their modeling. There are various bounds 
that can be studied, particularly in the context of the new generalized parton dynamics typical 
of the skewed distributions [HE]. All these relations can be analized in perturbation theory and 
studied using the Renormalization Group (RG), although a complete description of their pertur- 
bative dynamics is still missing. This study, we believe, may require considerable theoretical effort 
since it involves a global understanding both of the (older) forward (DGLAP) dynamics and of 
the generalized new dynamics encoded in the skewed distributions. 

In this context, a program aimed at the study of various bounds in perturbation theory using 
primarily a parton dynamics in x-space has been outlined [6j. This requires accurate algorithms 
to solve the equations up to next-to-leading order (NLO). Also, underlying this type of description 
is, in many cases, a probabilistic approach [3J which has some interesting consequences worth of 
a closer look . In fact, the DGLAP equation, viewed as a probabilistic process, can be rewritten 
in a master form which is at the root of some interesting formal developements. In particular, a 
wide set of results, available from the theory of stochastic processes, find their way in the study 
of the evolution. Two of us have elaborated on this issue in previous work and proposed a 
Kramers- Moyal expansion of the DGLAP equation as an alternative way to describe its dynamics. 
Here, this analysis will be extended to the case of the non-forward evolution. 

With these objectives in mind, in this study we test x-space algorithms up to NLO and verify 
their accuracy using a stringent test: Soffer's inequality. As usual, we are bound to work with spe- 
cific models of initial conditions. The implementations on which our analysis are based are general, 
with a varying flavour number n/ at any threshold of intermediate quark mass in the evolution. 
Here, we address Soffer's inequality using an approach based on the notion of "superdistributions" 
|3], which are constructs designed to have a simple (positive) evolution thanks to the existence 
of an underlying master form E|. The original motivation for using such a master form (also 
termed kinetic or probabilistic) to prove positivity has been presented in while further exten- 
sions of these arguments have been presented in [6J. In a final section we propose the extension of 
the evolution algorithm to the case of the skewed distributions, and illustrate its implementation 
in the non-singlet case. As for the forward case, numerical tests of the inequality are performed 
for two different models. We show that even starting from a saturated inequality at the lowest 
evolution scale, the various models differ significantly even for a moderate final factorization scale 
of Q = 100 GeV. Finally, we illustrate in another application the evolution of the tensor charge 
and show that, in the models considered, differences in the prediction of the tensor charge are 
large. 



2 Prelude to x-space: A Simple Proof of Positivity of hi to 
NLO 

There are some nice features of the parton dynamics, at least in the leading logarithmic approx- 
imation (LO), when viewed in x-space, once a suitable "master form" of the parton evolution 
equations is identified. 
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evolution equation itseii. ine topic lias ueen auuresseu oeiore in lv; j^j anu reanaiizeu in more 
detail in |J| where, starting from a kinetic interpretation of the evolution, a differential equation 
obtained from the Kramers-Moyal expansion of the DGLAP equation has also been proposed. 

The arguments of refs. |U E] are built around a form of the evolution equation which has a 
simple kinetic interpretation and is written in terms of transition probabilities constructed from 
the kernels. 

The strategy used, at least in leading order, to demonstrate the positivity of the LO evolution 
for special combinations of parton distributions Q ± j4j, to be defined below, or the NLO evolution 
for hi, which we are going to address, is based on some results of ref.|l], briefly reviewed here, in 
order to be self-contained. 

A master equation is typically given by 

d f 

f(x, t) = / dx (w(x\x )f(x, r) — w(x \x)f(x, r)) dx (1) 



Or 

and if through some manipulations, a DGLAP equation 



dq(x, Q 2 ) f 1 dy 2 

P(x/y)q{y,Q ), (2) 



d\og(Q 2 ) Jx y 
with kernels P(x), is rewritten in such a way to resemble eq. (QJ) 

U{x, r) = I' dyP (A^A_r d JLp (V) !M, (3 ) 
dr Jx \y J y Jo y \xJ x 

with a (positive) transition probability 

w{x\y) = ^P(x/y) V — ' (4) 

then positivity of the evolution is established. 

For equations of non-singlet type, such as those evolving = q — q, the valence quark 
distribution, or hi, the transverse spin distribution, this rewriting of the equation is possible, 
at least in LO. NLO proofs are, in general, impossible to construct by this method, since the 
kernels turn out, in many cases, to be negative. The only possible proof, in these cases, is just 
a numerical one, for suitable (positive) boundary conditions observed by the initial form of the 
parton distributions. Positivity of the evolution is then a result of an unobvious interplay between 
the various contributions to the kernels in various regions in x-space. 

In order to discuss the probabilistic version of the DGLAP equation it is convenient to separate 
the bulk contributions of the kernels (x < 1) from the edge point contributions at x — 1. For this 
purpose we recall that the structure of the kernels is, in general, given by 

P(z) = P{z) - 5(1 - z) f 1 P(z) dz, (5) 

Jo 

where the bulk contributions (z < 1) and the edge point contributions (~ S(z — 1)) have been 
explicitely separated. We focus on the transverse spin distributions as an example. With these 
prerequisites, proving the LO and NLO positivity of the transverse spin distributions is quite 
straightforward, but requires a numerical inspection of the transverse kernels. Since the evo- 
lutions for A T q^> = h\ are purely non-singlet, diagonality in flavour of the subtraction terms 
(~ Jq w(y\x)q(x, r)) is satisfied, while the edge-point subtractions can be tested to be positive 
numerically. We illustrate the explicit construction of the master equation for hi in LO, since 
extensions to NLO of this construction are rather straighforward. 
In this case the LO kernel is given by 



2 



A T P£\x) = C F 



xU 2 



2 + -6(l-x) 



(6) 



and by some simple manipulations we can rewrite the corresponding evolution equation in a 
suitable master form. That this is possible is an elementary fact since the subtraction terms can 
be written as integrals of a positive function. For instance, a possibility is to choose the transition 
probabilities 

wi[x\y] = ^L^-j--2 y je(y>x)e(y<l) 



w 2 [y\x} = ^ [ i-y/ x ~ I) 6 ( y > ~ x ^ y < 0) 

(7) 

which reproduce the evolution equation for hi in master form 

dhi f l f l 

— = dywi(x\y)hi(y,r) - dyw 2 (y\x)h 1 (x,r). (8) 
dr Jo Jo 

The NLO proof of positivity is also rather straightforward. For this purpose we have analized 

numerically the behaviour of the NLO kernels both in their bulk region and at the edge-point. 

We show in Table 1 of Appendix B results for the edge point contributions to NLO for both of 

the A.tP± components, which are numerically the same. There we have organized these terms 

in the form ~ C8{\ — x) with 

C=-log(l-A)A + B, (9) 

with A and B being numerical coefficients depending on the number of flavours included in the 
kernels. The (diverging) logarithmic contribution (~ Jq dz/(\—z)) have been regulated by a cutoff. 
This divergence in the convolution cancels when these terms are combined with the divergence 
at x = 1 of the first term of the master equation (JBJ) for all the relevant components containing 
"+" distributions. As for the bulk contributions (x < 1), positivity up to NLO of the transverse 
kernels is shown numerically in Fig. (0). All the conditions of positivity are therefore satisfied 
and therefore the Ar±q distributions evolve positively up to NLO. The existence of a master form 
of the equation is then guaranteed. 

Notice that the NLO positivity of Ax±q implies positivity of the nucleon tensor charge [S] 

5q = jT 1 dx (h\(x) - hf x) ) (10) 

for each separate flavour for positive initial conditions. As we have just shown, this proof of 
positivity is very short, as far as one can check numerically that both components of eq.(JHJ) are 
positive. 



3 Soffer's inequality 

Numerical tests of Soffer's inequality can be performed either in moment space or, as we are 
going to illustrate in the next section, directly in x-space, using suitable algorithms to capture the 
perturbative nature of the evolution. We recall that Soffer's inequality 

\hi(x)\<q + (x) (11) 

sets a bound on the transverse spin distribution hi(x) in terms of the components of the positive 
helicity component of the quarks, for a given flavour. An original proof of Soffer's inequality in 
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interpretation 01 tne evolution equations. 

We recall that h±, also denoted by the symbol 



A T q(x, Q 2 ) = q^(x, Q 2 ) - q l (x, Q 2 ), (12) 

has the property of being purely non-singlet and of appearing at leading twist. It is identifiable in 
transversely polarized hadron-hadron collisions and not in Deep Inelastic Scattering (from now on 
we will omit sometime the x-dependence in the kernels and in the distributions when obvious). In 
the following we will use interchangeably the notations h\ = h\ and Axq to denote the transverse 
asymmetries. We introduce also the combinations 

A T (q + q) = hl + hl 
A T q^=A T (q-q) = h\ - h\ 

A T q^ = £A T ( ft + fc) 

(13) 

where we sum over the flavor index (i), and we have introduced singlet and non-singlet contribu- 
tions for distributions of fixed helicities 

,(+) 



i 

Q+ ] = q+i - g+i = s. 



(14) 



In our analysis we solve all the equations in the helicity basis and reconstruct the various helicities 
after separating singlet and non-singlet sectors. We mention that the non-singlet sector is now 
given by a set of 2 equations, each involving ± helicities and the singlet sector is given by a 4-by-4 
matrix. 

In the singlet sector we have 

+Pf + ®G + + P q + G ®G-), 

dq( - ){x) -^(P + -^i +) +P ++ ^ +) 
+Pt G ®G + + Pf + ®G-), 

a ° ' ^Gq « „(+) , „G, « (+) 



dt 




dG + ( 


x) 


dt 




dG-( 


x) 


dt 



2n (PX®qr + P+*®Q'- 
+P GG ®G + + P GG ®G_), 

^{p?-® q { : ] +p? + ® q M 

+P GG ®G + + P GG ®G_). (15) 
while the non-singlet (valence) analogue of this equation is also easy to write down 

d -^r 1 = £ (p +- q + ] + p ++ (16) 
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tne nencity. 

Similarly to the unpolarized case the flavour reconstruction is done by adding two additional 
equations for each flavour in the helicity ± 

x ± ,i = <&> - (17) 



whose evolution is given by 

d£}(x) 



1 7 ■ 



dt 2tt 



dt 2n 



(18) 



The reconstruction of the various contributions in flavour space for the two helicities is finally 
done using the linear combinations 

q±i = I (q£ + X ±l + ^<?i +) ) ■ (19) 

We will be needing these equations below when we present a proof of positivity up to LO, and 
we will thereafter proceed with a NLO implementation of these and other evolution equations. 
For this we will be needing some more notations. 

We recall that the following relations are also true to all orders 

P(x) = ±(P ++ (x) + P + _(x)) 
= i(P-(x) + P_ + (x)) 

between polarized and unpolarized (P) kernels and 

P ++ (x)=P^(x), P_ + {x) = P + _{x) (20) 

relating unpolarized kernels to longitudinally polarized ones. Generically, the kernels of various 
type are expanded up to NLO as 

P(x) = ^\x) + {^-) 2 P^\x\ (21) 
and specifically, in the transverse case we have 

A T P« ± = A T P«±A T P«, (22) 

(23) 

with the corresponding evolution equations 

d 



dlnQ- 



AtQ±(Q ) = A T P qq , ± (a s (Q 2 )) ® A T q±(Q 2 ) . (24) 
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p(0) 
^NS±,+ 

p(0) 
p(0) 



P£] ++ = P' v >=N c [—^— + --l-x-x 2 )+f3 5(l-x) 



P 9 ( 9 0) +- = N c [3x + --3-x 2 \ . (25) 




An inequality, such as Soffer's inequality, can be stated as positivity condition for suitable 
linear combinations of parton distributions |4| and this condition can be analized - as we have just 
shown for the hi case - in a most direct way using the master form. 

For this purpose consider the linear valence combinations 

Q+ = Q+ + h 
Q = q+ - h 

(26) 

which are termed "superdistributions" in ref.|l]. Notice that a proof of positivity of the Q dis- 
tributions is equivalent to verify Soffer's inequality. However, given the mixing of singlet and 
non-singlet sectors, the analysis of the master form is, in this case, more complex. As we have 
just mentioned, what can spoil the proof of positivity, in general, is the negativity of the kernels 
to higher order. We anticipate here the result that we will illustrate below where we show that a 
LO proof of the positivity of the evolution for Q can be established using kinetic arguments, being 
the kernels are positive at this order. However we find that the NLO kernels do not satisfy this 
condition. In any case, let's see how the identification of such master form proceeds in general. 
We find useful to illustrate the result using the separation between singlet and non-singlet sectors. 
In this case we introduce the combinations 

Qi +) = qi +) ±h[ +) 

(27) 

with h [ i ] = A T ? (±) . 

Differentiating these two linear combinations (|2Tjl we get 



rfio g (g 2 ) 



(28) 



which can be rewritten as 



dio g (g 2 ) 

d Q+ } _ 1 (jd(~) o(-)\nH i 1 /d(-) , nH , d(-)„(-) 

dlog(Q 2 ) 



(29) 
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Figure 1: Plot of the transverse kernels. 

with P(~) = P NS being the non-singlet (NS) kernel. 
At this point we define the linear combinations 



and rewrite the equations above as 

dQ + i 



dhg(Q 2 ) 

dQi+ 
rflog(Q 2 ) 



P?± = 2 ( P ++ ± p r) 



(30) 



(31) 



where we have reintroduced i as a flavour index. From this form of the equations it is easy to 
establish the leading order positivity of the evolution, after checking the positivity of the kernel 
and the existence of a master form. 

The second non-singlet sector is defined via the variables 



Xi± - Qi± Qi± 

which evolve as non-singlets and the two additional distributions 

Qx«,± = Xi+ ± h % i +) . 
Also in this case we introduce the kernels 

P?l = \(P ++ ±A T PM) 

to obtain the evolutions 

dQ X i + 



(32) 
(33) 
(34) 



rflog(Q 2 ) 

dQ i+ 
dlog(Q 2 ) 



P++Q\i+ + P+-Q\i- + P+-Xi 



^X+-Qx»+ + P++Q\i- + P+-X 



(35) 
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to tne singiet equation 01 trie nencity oasis, using tne equations auove, tne uistnoutions ^± are 
then reconstructed as 

Q*± = ~ (q£ + Q&i + (36) 

and result positive for any flavour if the addends are positive as well. However, as we have just 
mentioned, positivity of all the kernels introduced above is easy to check numerically to LO, 
together with their diagonality in flavour which guarantees the existence of a master form. 

As an example, consider the LO evolution of Q±. The proof of positivity is a simple con- 
sequence of the structure of eq. (pTT]) . In fact the edge-point contributions appear only in 
i.e. they are diagonal in the evolution of Q±. The inhomogenous terms on the right hand side of 
(|3TJ) . proportional to g_ are are harmless, since the P + _ kernel has no edge-point contributions. 
Therefore under 1) diagonality in flavour of the subtraction terms and 2) positivity of first and 
second term (transition probabilities) we can have positivity of the evolution. A refined arguments 
to support this claim has been presented in [H]. 

This construction is not valid to NLO. In fact, while the features of flavour diagonality of the 
master equation are satisfied, the transition probabilities w(x,y) are not positive in the whole 
x, y range. The existence of a crossing from positive to negative values in P+ + can, in fact, be 
established quite easily using a numerical analysis. We illustrate in Figs. J2J and (JHJ) plots of the Q 
kernels at LO and NLO, showing that, at NLO, the requirement of positivity of some components 
is violated. The limitations of this sort of proofs -based on kinetic arguments- are strictly linked 
to the positivity of the transition probabilities once a master form of the equation is identified. 



4 An x-space Expansion 

We have seen that NLO proofs of positivity, can be -at least partially- obtained only for suitable 
sets of boundary conditions. To this purpose, we choose to investigate the numerical behaviour 
of the solution using x-space based algorithms which need to be tested up to NLO. 

Our study validates a method which can be used to solve evolution equations with accuracy in 
leading and in next-to-leading order. The method is entirely based on an expansion jTH] used in 
the context of spin physics |H] and in supersymmetry An interesting feature of the expansion, 
once combined with Soffer's inequality, is to generate an infinite set of relations among the scale 
invariant coefficients (A nj B n ) which characterize it. 

In this approach, the NLO expansion of the distributions in the DGLAP equation is generically 
given by 

where, to simplify the notation, we assume a short-hand matrix notation for all the convolution 
products. Therefore f(x, Q 2 ) stands for a vector having as components any of the helicities of the 
various flavours (Q±,q±,G±). The ansatz implies a tower of recursion relations once the running 
coupling is kept into account and implies that 

A n+1 (x) = -J-P (0) ® A n (x) (38) 



Po 



to leading order and 

B n+ i{x) = -B n (x) 
= -BJx 



^- 2 P^®A n {x) 
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(39) 

to NLO, relations which are solved with the initial condition B (x) = 0. The initial conditions 
for the coefficients A (x) and B (x) are specified with q(x, Ql) as a leading order ansatz for the 
initial distribution 

A (x) = 5(l-x)®q(x,Ql) = q (x) (40) 
which also requires B (x) = 0, since we have to satisfy the boundary condition 

Aq{x) + a B (x) = q (x). (41) 

Again, other boundary choices are possible for A (x) and B (x) as far as l|ITl) is fullfilled. 
If we introduce Rossi's expansion for hi, q + , and the linear combinations Q± (in short form) 



h ~ (A h n ,B h n 
q ± ~ {A*,B?\ 

(42) 

we easily get the inequalities 

(-ir(^++4)>0 (43) 

and 

(A? - A*) >0 (44) 

valid to leading order, which we can check numerically. Notice that the signature factor has to be 
included due to the alternation in sign of the expansion. To next to leading order we obtain 



{-l) n+1 (Al+{x)+a{Q 2 )B q n + {x)) <(-ir(A h n (x) + a(Q 2 ) B h n (x)) < (-1)" (A*+{z) + a(Q 2 )B^(x)) , 

(45) 

valid for n > 1, obtained after identification of the corresponding logarithmic powers log (a(Q 2 )) 
at any Q. In general, one can assume a saturation of the inequality at the initial evolution scale 

Q_(x, Ql) = h^x, Ql) - ~q+(x, Ql) = 0. (46) 

This initial condition has been evolved in Q solving the equations for the Q± distributions to 
NLO. 



5 Relations among moments 

In this section we elaborate on the relation between the coefficients of the recursive expansion as 
defined above and the standard solution of the evolution equations in the space of Mellin moments. 
We will show that the two solutions can be related in an unobvious way. 

Of our concern here is the relation between the Mellin moments of the coefficients appearing 
in the expansion 

A(N) = C dxx N - 1 A(x) 
Jo 

B(N) = [ dxx N - 1 B(x) 
Jo 

(47) 
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A T q±(N,Q 2 ) = f 1 dxx N ~ l A T (x,Q 2 )). (48) 
Jo 

For this purpose we recall that the general (non-singlet) solution to NLO for the latter moments 
is given by 

A T g ± (iV,Q 2 ) = K(QlQ\N)(^/) A T q ± (N,Q 2 ) 
with the input distributions A T q±(Ql) at the input scale Q and where we have set 



(49) 



In the expressions above we have introduced the corresponding moments for the LO and NLO 
kernels (A T P$> N , A T P«f). 

We can easily get the relation between the moments of the coefficients of the non-singlet x-space 
expansion and those of the parton distributions at any Q, as expressed by eq. (jlHjl 

A n (N) + a s B n (N) = A T q ± (N, Q 2 )K(Q , Q, N) ( Z ^^ i ^) ■ (50) 

As a check of this expression, notice that the initial condition is easily obtained from (IHOI) 
setting Q — > Qo, n ~ > 0, thereby obtaining 

A» 8 {N) + a s B^ s (N) = A T q ± (N, Q 2 ) (51) 

which can be solved with A$ S (N) = A T q±(N, Q 2 ) and Bg s (N) = 0. 

It is then evident that the expansion lpT7j) involves a resummation of the logarithmic contribu- 
tions, as shown in eq. ljH0~J) . 

In the singlet sector we can work out a similar relation both to LO 



with 



ei = (P {0 \N) - \ 2 1 

A\ — A2 v 

e 2 = — ^-(-P<°>(JV)+Ail 

A 2 — Ai v 



Al,2 

and to NLO 
where 



I ( P (o) (AT) + p(J (AT) ± y (pW(iV) _ pW(iV)) 2 + 4P£\N)P£\N) ) , (53) 



A„(iV) + a sJ B n (iV) = Xif-^ i ) +X 2 f-^) (54) 



a /-2 e 2 Rei 
Xi = ei + — — eiRei + 



2vr V/?o Ai — A 2 — Po/2, 

a (-2 t eiRe 2 

X2 = e 2 + — — e 2 Re 2 



2vr V/5o Aa-Ai-Aj/2, 



(55) 
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Figure 2: Plot of the LO kernels for the Q distributions 



with 



K = p( 1 \N) - Ap(°)(/y). 

2po 



(56) 



Notice that A n (N) and B n (N), P (0) (iV), pW(N), in this case, are all 2-by-2 singlet matrices. 
Prior to discuss some phenomenological application of this method, left to a final section, we will 
now illustrate its application to the case of the evolution of hi up to NLO. Its extension to the 
case of the generalized parton distributions will follow afterwards. 



6 An Example: The Evolution of the Transverse Spin Dis- 
tributions 

LO and NLO recursion relations for the coefficients of the expansion can be worked out quite easily, 
although the numerical implementation of these equations is far from being obvious. Things are 
somehow simpler to illustrate in the case of simple non-singlet evolutions, such as those involving 
transverse spin distributions, as we are going to show below. Some details and definitions can be 
found in the appendix. 

For the first recursion relation (eq. (JHBj!) we have 



'/V 



--A T P$(x) ® A±(x) 
dyyA±(y) -xA±(x) 



x y 



y-x 



+ A±(x) log(l 



x) 



+ 



<)(jfTH +c '(-5)I**> 



(57) 



As we move to NLO, it is convenient to summarize the structure of the transverse kernel ArP^ 1 ' 
as 



A T P±,(1) 



x) 



Kf(x)5{\ -x) + K±(x)S 2 {x) + Kf(x) log(x) 
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Hence, for the (+) case we have 



x y 



K+(z)S 2 (z)+K+(z)\og(z) 



+ \og\z)Ki{z) + log(*) log(l - z)Kt{z) A+{y) + 




dyyA+(y) -xA+(x) 



x y 



y-x 



+ AZ(x) log(l -x)\+Kt I —A+{y) 

x y 



dy 



(59) 



where z = x/y. For the (— ) case we get a similar expression. 

Now we are ready to write down the expression for the B n+1 (x) coefficient to NLO, similarly 
to eq. (|39|) . So we get (for the (+) case) 



-B+(x) + 



2PZ 



2C f 



dyyA+(y) - xA+(x 



-2C F U^A+{y)\+C F \A+{x) 




KtA + n (x 



+ K+(z) \og(z) + log 2 (z)K+(z) + log(z) log(l - z)K+{z) 




+ 



GA- 




dyyA+(y) -xA+(x) 



x y 



y-x 



+ A+(x) log(l- 



x , 



Jx y 



1 dyyB±(y) - xB±(x) 



Bt(x) logfl -x] 



x y 



y-x 



+ 



-S)(/:?H + -B)^- 



(60) 



As we have already mentioned, the implementation of these recursion relations require particular 
numerical care, since, as n increases, numerical instabilities tend to add up unless high accuracy is 
used in the computation of the integrals. In particular we use finite element expansions to extract 
analitically the logarithms in the convolution (see the discussion in Appendix A). NLO plots of 
the coefficients A n (x) + a(Q 2 )B n (x) are shown in figs. (I4|5j) for a specific set of initial conditions 
(GRSV, as discussed below). As the index n increases, the number of nodes also increases. A 
stable implementation can be reached for several thousands of grid-points and up to n ~ 10. Notice 
that the asymptotic expansion is suppressed by n! and that additional contributions (n > 10) are 
insignificant even at large (> 200 GeV) final evolution scales Q. 



7 Nonforward Extensions 

In this section we finally discuss the nonforward extension of the evolution algorithm. In the case of 
nonforward distributions a second scaling parameter ( controls the asymmetry between the initial 
and the final nucleon momentum in the deeply virtual limit of nucleon Compton scattering. The 
solution of the evolution equations, in this case, are known in operatorial form. Single and double 
parton distributions are obtained sandwiching the operatorial solution with 4 possible types of 
initial/final states < >,< j?| — 1 >, < p'\...\p >, corresponding, respectively, to the case of 

diagonal parton distributions, distribution amplitudes and, in the latter case, skewed and double 
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Figure 3: Plot of the NLO kernels for the Q distributions, showing a negative behaviour at large 
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parton distributions f2]. Here we will simply analize the non-singlet case and discuss the extension 
of the forward algorithm to this more general case. Therefore, given the off-forward distributions 
H q (x,£), in Ji's notation, we set up the expansion 



k=0 



^ log *(^W)£^ log *(4« 2 

k\ \a(Qo) f n k\ \a(^ 0) 



(61) 



which is the natural extension of the forward algorithm discussed in the previous sections. We 
recall that in the light-cone gauge H(x,£) is defined as 



H q (x,^A 2 )) = - 



1 rdy- 



2tt 



- p+ y- (p'imo, y -, o ± )J T +w(o, y -, o ± )\p) 



(62) 



with A = P' - P, P+ = 1/2(P + P) H (symmetric choice) and £P = 1/2 A+. 

This distribution describes for x > £ and x < — £ the DGLAP-type region for the quark and the 
antiquark distribution respectively, and the ERBL [22J (see also [23J for an overview) distribution 
amplitude for — £ < x < £. In the following we will omit the A dependence from H q . 

Again, once we insert the ansatz (|6Tf into the evolution equations we obtain an infinite set 
of recursion relations which we can solve numerically. In LO, it is rather simple to relate the 
Gegenbauer moments of the skewed distributions and those of the generalized scaling coefficients 
A n . We recall that in the nonforward evolution, the multiplicatively renormalizable operators 
appearing in the light cone expansion are given in terms of Gegenbauer polynomials [2J. The 
Gegenbauer moments of the coefficients A n of our expansion (|6T|) can be easily related to those of 
the off-forward distribution 



C n (£, Q 2 ) = C £ C 3 J 2 (z/0H(z, e, Q 2 )di 
The evolution of these moments is rather simple 

->2\ \ 7n/A) 



C n (CQ 2 )=C n (CQ 2 ) 



a(Q 2 

AQl) 



(63) 



(64) 
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7n 



(n + l)(n + 2) J 



(65) 



being the non-singlet anomalous dimensions. If we define the Gegenbauer moments of our expan- 



sion 



4%, Q 2 ) = C £ CW(z/t)H(zt, Q 2 )dz 
we can relate the moments of the two expansions as 



7n 



(66) 



(67) 



Notice that expansions similar to (|HT1) hold also for other choices of kinematical variables, such 
as those defining the non- forward distributions j2j, where the t-channel longitudinal momentum 
exchange A + is related to the longitudinal momentum of the incoming nucleon as A = (P. We 
recall that H q (x.£) as defined in [T] can be mapped into two independent distributions F q (X, () 
and T q {X, () through the mappings [E| 



x x 
x 2 

^(x 2 ,0 



(i + 

(i + 
C/(2-C) 

-H q (x 1 ,£) 



C/2 
-1 



C/2 



(68) 



in which the interval — 1 < x < 1 is split into two coverings, partially overlapping (for — £ < x < £, 
or ERBL region) in terms of the two variables — £ < x\ < 1 (0 < X\ < 1) and — 1 < x 2 < £ 
(0 < X2 < 1). In this new parameterization, the momentum fraction carried by the emitted quark 
is X, as in the case of ordinary distributions, where it is parametrized by Bjorken x. For definitess, 
we focus here on the DGLAP-like (X > () region of the non-singlet evolution. The non-singlet 
kernel is given in this case by (x = X) 



*dx,C) = -C F ( 

7r \y — x 



XX 

1 + — 

yy 



-5{x-y) [ dz \ + Z ) 
Jo I — z J 



we introduce a LO ansatz 



k=0 



k\ 




(69) 



(70) 



and insert it into the evolution of this region to obtain the very simple recursion relations 



A n+1 (X,C) 



2 

To 



dyyA n (y,()-xA n (X,() 2 



x y 

c F Mx,0 



+ In 



-X 

[l-Xf(l 



Po 



dy(X-Q (yA n (X,()-XA n (y,Q) 

* y(y-Q y~ x 



i-C 



(71) 



The recursion relations can be easily reduced to a weighted sum of contributions in which ( is 
a spectator parameter. Here we will not make a complete implementation, but we will illustrate 
in an appendix the general strategy to be followed. There we show a very accurate analytical 
method to evaluate the logarithms generated by the expansion without having to rely on brute- 
force computations. 
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Positivity of the non-singlet evolution is a simple consequence of the master-form associated to 
the non-forward kernel (JHSJl- As we have already emphasized above, positivity of the initial 
conditions are sufficient to guarantee a positivity of the solution at any scale Q. The master- form 
of the equation allows to reinterpret the parton dynamics as a random walk biased toward small-x 
values as r = log(Q 2 ) increases. 

In the non-forward case the identification of a transition probability for the random walk [H] 
associated with the evolution of the parton distribution is obtained via the non-forward transition 
probability 



W((x\y) 
w' c (y\x) 



a 1 

— Cyj? 

7T y — x 

a: 



1 + 



x {x - C) 

y y-( 



6{y > x) 



Co t + y \ o(v < x) 



it x 2 (x — y) 
and the corresponding master equation is given by 



(72) 



dT q 
Ik 



dyw^x^FqiyX; 



T 



dywUy\x)F q (x,(,r), 



(73) 



that can be re-expressed in a form which is a simple generalization of the formula for the forward 
evolution [0] 



dT q <> 



d log Q 2 Jx 



dyw c (x\y)T q (y,(^)- / dy w[{y\x)T q {x, (, r) 

<j o 

a(x) r—x 



ra(x) r—x 

= -J dyw < (x + y\x)*J 7 q (x,C,r) + j dyw' c (x + y\x)J 7 q (x,C,r), (74) 
where a Moyal-like product appears 

—y( d x + x \ 

w^(x + y\x) * T q {x, (, t) = w^(x + y\x)e \ ' Fqipz-, (75) 

and a(x) = x — 1. A Kramers-Moyal expansion of the equation allows to generate a differential 
equation of infinite order with a parametric dependence on ( 

dT q r° 



We define 



dyw c (x + y\x)J 7 q (x,(,r)+ dyw' c (x + y\x)F q (x,(,r) 

diOgQ 1 Ja(x) JO S 

~ ra(x) {-y) n 

E / dy K -^-d x n (w c (x + y\x)F q (x, C, r)) . (76) 

n=l J0 n - 



rO r — x 

d (x,C) = / dyw i (x + y\x)J 7 q (x,(,T)+ dyw'Ax + y\x)F q {x, C, r) 

J a(x) JO 
ra(x) 

a n (x,C) = / dyy n wc:(x + y\x)T q {x, £t) 

J 

a n (x,() = dyy n d x n (w c (x + y\x)F q (x,C,T)) n = 1,2,... (77) 

If we arrest the expansion at the first two terms (n = 1, 2) we are able to derive an approximate 
equation describing the dynamics of partons for non-diagonal transitions. The procedure is a 
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purpose we use trie identities 



ai(x,() = d x ai(x, () — a(x)d x a(x)w^(x + a(x)\x)J 7 q (x, (, r) 
02(^,0 = dla 2 (x, C) — 2a(x)(d x a(x)) 2 w^(x + a(x)\x)J 7 q (x, (, r) 
— a(a;) 2 9 a: Q;(a;)5 a ; (u^(a; + t)) 

-a 2 (x)9 x a(a;)9 ;c (iu ( (a; + y|a;)J^(a;, C, r)) |y= a (x)- (78) 

which allow to compute the first few coefficients of the expansion. Using these relations, the 
Fokker-Planck approximation to this equation can be worked out explicitely. We omit details 
on the derivation which is unobvious since particular care is needed to regulate the (canceling) 
divergences and just quote the result. 
A lengthy computation gives 

dT q a ( S0-3 xq-i \ . . 

-*F = ^HmF + F0 + x T 5(x,c ' t) 

, a n ( Xl -3 ^i-i \ a t ( /• \ , a n x 0-3 n2 T / r \ 

+~Cf (^7—^ + ) d x F q (x, C, r) + -C F ^-^^(x, C , r) 



(79) 



where we have defined 

- ((-1 + xf (17a; 3 - C 2 (3 + 4C) + 3x( (3 + 5C) - 3a; 2 (3 + 70)) 



XO-3 



12 x s 

-29a; 4 - 3 + x 2 (-1 + + 2C - 2a; (1 + 30 + a; 3 (12 + 230 



3a; 3 



1 3 rtl (1-x) 
^0,0 = 4 + — -- + 21og v 



2a; 2 x x 

- ((-1 + 6a; - 15a; 2 + 14a; 3 ) (x - ()) 

3x 2 

1 5a; r o 23a; 4 7( 3( 5a< 

a;i _q = h 5a; 1 1 

' 2 3 6 3 Ax 2 

1r 2A 13la; 3 C 5C 2 ( 2 10 , 2 39x 2 ( 2 OA , 8a< 3 

- f(-l + a;) 2 (a; - C) 2 (3 + 23a; 2 + 4( - 2x {7 + 80)) 

^2,-3 = ^- (80) 

This equation and all the equations obtained by arresting the Kramers-Moyal expansion to 
higher order provide a complementary description of the non-forward dynamics in the DGLAP 
region, at least in the non-singlet case. Moving to higher order is straightforward although the 
results are slightly lengthier. A full-fledged study of these equations is under way and we expect 
that the DGLAP dynamics is reobtained - directly from these equations - as the order of the 
approximation increases. 



9 Model Comparisons, Saturation and the Tensor Charge 

In this last section we discuss some implementations of our methods to the standard (forward) 
evolution by doing a NLO model comparisons both in the analysis of Soffer's inequality and for 
the evolution of the tensor charge. We have selected two models, motivated quite independently 
and we have compared the predicted evolution of the Soffer bound at an accessable final evolution 
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tins point we recan tnat 111 order to generate suitauie initial conditions tor tne analysis 01 ooner s 
inequality, one needs an ansatz in order to quantify the difference between its left-hand side and 
right-hand side at its initial value. 

The well known strategy to build reasonable initial conditions for the transverse spin distribu- 
tion consists in generating polarized distributions (starting from the unpolarized ones) and then 
saturate the inequality at some lowest scale, which is the approach we have followed for all the 
models that we have implemented. 

Following Ref. [IB] (GRSV model), we have used as input distributions - in the unpolarized 
case - the formulas in Ref. [P7j . calculated to NLO in the MS scheme at a scale Ql = 0.40 GeV 2 

x(u-u)(x,Ql) = 0.632x°' 43 (l -x) 3 09 (l + 18.2x) 

x{d-d){x,Q 2 ) = 0.624(1 - x) ho x(u -u)(x,Ql) 

x(d-u)(x,Q 2 ) = 0.20a; - 43 (l -x) 12 - 4 (l - 13.3v^ + 60.0a;) 

x(u + d)(x,Ql) = 1.24x a20 (l -x) 85 (l -2.3v^ + 5.7x) 

xg{x,Ql) = 20.80x L6 (l - x) 4A (81) 

and xqi(x, Ql) = xqi(x, Ql) = for qi = s, c, b, t. 

We have then related the unpolarized input distribution to the longitudinally polarized ones 
by as in Ref. |THj 

xAu{x,Ql) = 1.019x - 52 {l-x)°- 12 xu{x,Ql) 
xAd(x,Q 2 ) = -0.669x 0A3 xd(x,Q 2 ) 
xAu(x,Ql) = -0.272a: - 38 xu(x, Q 2 ) 
xAd(x, Ql) = xAu(x, Qq) 

xAg(x,Q 2 ) = 1.419x L43 (l - x) 0A5 xg(x,Ql) (82) 

and xAqi(x, Ql) = xAqj{x, Ql) = for q^ = s, c, b, t. 

Following [21], we assume the saturation of Soffer inequality: 

a ( /o2\ x Qi( x i Ql) + xAqi(x, Ql) 
xA T qi(x,Q ) = (83) 

and study the impact of the different evolutions on both sides of Soffer's inequality at various final 
evolution scales Q. 

In the implementation of the second model (GGR model) we have used as input distributions 
in the unpolarized case the CTEQ4 parametrization [11], calculated to NLO in the MS scheme at 
a scale Qo = 1-0 GeV 



x(u-u){x,Ql) = 1.344x a501 (l-x) 3 - 689 (l + 6.402x a873 ) 

x(d-d)(x,Q 2 ) = 0.64x - 501 (l-x) 4 - 247 (l + 2.69a; a333 ) 

xs(x, Ql) = xs(x, Ql) = 0.064x- ai43 (l - x) 8 - 041 (l + 6.112a;) 

x(d-u)(x,Ql) = 0.071a; a501 (l-x) 8 - 041 (l + 30.0x) 

x(u + d)(x,Ql) = 0.255a;- - 143 (l-x) 8 041 (l + 6.112x) 

xg(x,Q 2 ) = 1.123x-°- 206 (l-a;) 4 - 673 (l + 4.269x L508 ) (84) 

and xqi(x,Ql) = xqi(x,Ql) = for q^ = c,b,t and we have related the unpolarized input distri- 
bution to the longitudinally polarized ones by the relations (2^ 
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xAu(x, Ql) 



xAu(x, Q ) 
xAd(x,Q 2 ) 
xAd(x, Ql) 
xAs(x, Ql) = xAs(x, Ql) 



xt] u (x)xu(x,QI) 
cos 9d(x,Q 
xr) d (x)xd(x, Qq 



cos 6 D (x, Ql 



x{u — u 
2 ) 

1 



-x(d — d) 



(x, Ql) + xAu(x, Q 



x(d — d) (x, Qqj 



+ xAd(x,Q 2 ) 



xrj s (x)xs(x,Ql) 



(85) 



and xAqi(x, Ql) = xAqi(x, Ql) = for = c, b, t. 



A so-called "spin dilution factor" as defined in 
given by 

2a s (Q 2 ) (1 



which appears in the equations above is 



cos Ql 



1 + 



x) 



X 



(86) 



In this second (GGR) model, in regard to the initial conditions for the gluons, we have made use 
of two different options, characterized by a parameter r\ dependent on the corresponding option. 
The first option, that we will denote by GGR1, assumes that gluons are moderately polarized 



xAg(x,Q 2 ) 
Vu(x) = 7} d (x) 
r) s (x) 



x ■ xg(x,Ql) 
-2.49 + 2.8^ 
-1.67 + 2.1^, 



(87) 



while the second option (GGR2) assumes that gluons are not polarized 

xAg(x,Ql) = 
Vu(x) = rjd(x) = -3.03 + 3.0\/x 
Vs ( x ) = -2.71 + 2.9-v/s. 

We have plotted both ratios A T / f + and differences (xf + — xA T f) for various flavours as a function 
of x. For the up quark, while the two models GGR1 and GGR2 are practically overlapping, the 
difference between the GGR and the GRSV models in the the ratio Atu/u + is only slightly 
remarked in the intermediate x region (0.1 — 0.5). In any case, it is just at the few percent level 
(Fig. ©), while the inequality is satisfied with a ratio between the plus helicity distribution and 
transverse around 10 percent from the saturation value, and above. There is a wider gap in the 
inequality at small x, region characterized by larger transverse distribution, with values up to 40 
percent from saturation. A similar trend is noticed for the x-behaviour of the inequality in the 
case of the down quark (Fig. In this latter case the GGR and the GRSV model show a more 
remarked difference, especially for intermediate x-values. An interesting features appears in the 
corresponding plot for the strange quark (Fig.(JBJ), showing a much wider gap in the inequality 
(50 percent and higher) compared to the other quarks. Here we have plotted results for the two 
GGR models (GGR1 and GGR2). Differently from the case of the other quarks, in this case 
we observe a wider gap between lhs and rhs at larger x values, increasing as x — > 1. In figs, 
©and (fT0|) we plot the differences (xf + — xA T f) for strange and charm and for bottom and top 
quarks respectively, which show a much more reduced evolution from the saturation value up to 
the final corresponding evolving scales (100 and 200 GeV). As a final application we finally discuss 
the behaviour of the tensor charge of the up quark for the two models as a function of the final 
evolution scale Q. We recall that like the isoscalar and the isovector axial vector charges defined 
from the forward matrix element of the nucleon, the nucleon tensor charge is defined from the 
matrix element of the tensor current 



(PS T \^ l5 X a iP\P,S T ) 



25q a (Ql) 



(P»S V T - P v StP) 



(89) 
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dt is me transverse spin. 

In fig. (fTTf we plot the evolution of the tensor charge for the models we have taken in exam. 
At the lowest evolution scales the charge is, in these models, above 1 and decreases slightly 
as the factorization scale Q increases. We have performed an evolution up to 200 GeV as an 
illustration of this behaviour. There are substantial differences between these models, as one can 
easily observe, which are around 20 percent. From the analysis of these differences at various 
factorization scales we can connect low energy dynamics to observables at higher energy, thereby 
distinguishing between the various models. Inclusion of the correct evolution, up to subleading 
order is, in general, essential. 

10 Conclusions 

We have illustrated the use of x-space based algorithms for the solution of evolution equations in 
the leading and in the next-to-leading approximation and we have provided some applications of 
the method both in the analysis of Soffer's inequality and in the investigation of other relations, 
such as the evolution of the proton tensor charge, for various models. The evolution has been 
implemented using a suitable base, relevant for an analysis of positivity in LO, using kinetic 
arguments. The same kinetic argument has been used to prove the positivity of the evolution of h\ 
and of the tensor charge up to NLO. In our implementations we have completely relied on recursion 
relations without any reference to Mellin moments. We have provided several illustrations of the 
recursive algorithm and extended it to the non- forward evolution up to NLO. Building on previous 
work for the forward evolution, we have presented a master-form of the non-singlet evolution of 
the skewed distributions, a simple proof of positivity and a related Kramers Moyal expansion, 
valid in the DGLAP region of the skewed evolution for any value of the asymmetry parameter 
(. We hope to return with a complete study of the nonforward evolution and related issues not 
discussed here in the near future. 
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11 Appendix A. Weighted Sums 

In this appendix we briefly illustrate the reduction of recursion relations to analytic expressions 
based on finite element decompositions of the corresponding integrals. The method allows to write 
in analytic forms the most dangerous integrals thereby eliminating possible sources of instabilities 
in the implementation of the recursion relations. The method uses a linear interpolation formula 
for the coefficients A n , B n which, in principle can also be extended to higher (quadratic) order. 
However, enugh accuracy can be achieved by increasing the grid points in the discretization. Notice 
that using this method we can reach any accuracy since we have closed formulas for the integrals. 
In practice these and similar equations are introduced analytically as functions in the numerical 
integration procedures. 

Below, we will use a simplified notation (X = x for simplicity). 

We define P(x, () = xP(x, () and A(x, () = xA(x) and the convolution products 




(90) 
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oruereu irom leit to rignt ^xo,xi,X2, x n , x n +ij witn xq = x anu x n +i = i ueing trie upper euge 
of the integration region. One constructs a rescaled array (x, x/x n , x/x 2 , x/xi, 1). We define 
Sj = x/xj, and s n+ i = x < s n < s n _i < ...Si < s — 1- We get 

' r x i+i dy ( x s 



i=0 7 ' T 



X 



. p -CM(y,C) 



At this point we introduce the linear interpolation 



%v,Q= i- 



A(xi,C) + 



Xi+1 Xi 



■A,C(x i+1 ) 



(91) 



(92) 



and perform the integration on each subinterval with a change of variable y— > x/y and replace 
the integral J(x,Q with its discrete approximation J/v(x) to get 



J N (x,Q = A(x )-^— f l< ^P{y,Q{y-s,) 

1 — Si J si V 



i=l 
N 



si J Sl y 

+ EA( Xi ,o—^— T -P(y)(y - s i+1 ) 



E%0^- r _1 -P(y,C)(y-^-i) 



(93) 

with the condition A(xn+i,() = 0- Introducing the coefficients W(x,x,() an d W(xi,x,(), the 
integral is cast in the form 



where 



J N (x, C) = W(x, x, C)A(x, O + E ^ O^i, 

1=1 



W(x,x,C) = l X —{v-8 X )P(y,Q, 
i — s\ Jsi y 

W(Xi,x,Q = 



(94) 



■Si — Sj+i Js 4+ i y 



Sj-i — Sj y 



dy 



(y-Si +1 )P(y,0 
(y- Sl ^)P(yX). 



(95) 



For instance, after some manipulations we get 



1 ^y vMyX) ^MX,() = Ino(x)AM+ £ ( j n . (;r) - Jnt,(x)) A,(xi)-ln(l-x)A,(x,C) 
2/ A i+i 



where 
I (x) 

Ji(x) 
3[(x) 
3i(x) 



(96) 



1 



1 



log(si) +log(l - Si) 
1 - S i+ i 



Si — Si+i 
1 

Si—1 Sj 
1 



log 
log 



1 - Si 
I - Sj 
. 1 - Si-! 



s i+ i log 



Si-i log 



1 S 7 ; S? 



i *i+l 



1 — Sj 



Si-1. 



+ Si-i 



1 - gi-1 

1 - Si 



, i = 2,3, ..JV 



1-si 



logsi. 



(97) 
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12 Appendix B 



n f 


A 


B 


3 


12.5302 


12.1739 


4 


10.9569 


10.6924 


5 


9.3836 


9.2109 


6 


7.8103 


7.7249 



Table 1. Coefficients A and B for various flavour, to NLO for ArP ggj ±. 



13 APPENDIX C 

Here we define some notations in regard to the recursion relations used for the NLO evolution of 
the transverse spin distributions. 

For the (+) case we have these expressions 



K+{x) = —C F {-2n f {3 + 4tt 2 ) + N c (51 + 44vr 2 - 216C(3)) + 9C F (3 - 4tt 2 + 48C(3)) 
K t{x) = ^ F (-2C F + N C)X (9g) 

Ktlx) = ^ (100) 
/tf (s) = ^ (101) 
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pa 

a 
+ 
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-0.2 



ii—i 

n=2 
n=3 
n=4 



-0.3 
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x 

Figure 5: Coefficients A n (x) + a s (Q 2 )B n , with n — 0, . . . , 4 for a final scale Q = 100 GeV for the 
quark down. 



K+(x) = --C F (10n f + iV c (-67 + 3tt 2 )) 

K t( x ) = lc F (10n f + N c (-67 + 3tt 2 )) 
9 



and for the (— ) case we have 
1 



(102) 

(103) 
(104) 



72 

2C F (+2C F - N c )x 
1 + x 

C F (9C F -llN c + 2n f )x 



Ki(x) = —C F (-2n f (3 + 4tt 2 ) + N c (51 + 44tt 2 - 216C(3)) + 9C F (3 - Air 2 + 48((3)) 
K;(x) 
K~(x) 

K-{x) 



3{x- 1) 



C F N c x 



l-x 

F 5 



AClx 



l-x 

K^{x) = -\c F {lQn } + iV c (-67 + 3vr 2 )) 

y 

iff (x) = ~C F (10n f - 18C F (x - 1) + N c (-76 + 3tt 2 + 9x)) 



(105) 
(106) 
(107) 
(108) 
(109) 



(110) 
(111) 
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Figure 7: Test of Soffer's inequality for quark down at Q = 100 GeV for different models 
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Figure 8: Test od Soffer's inequality for quark strange at Q = 100 GeV for different models 
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Figure 9: Soffer's inequality for strange and charm in the GRSV model. 
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Figure 10: Soffer's inequality for bottom and top in the GRSV model. 



25 



+ 




up(GRSV) — 
up(GGR) — ■ 

down(GRSV) 

down (GGR) 




















VVVVJVVVVVl 



20 40 60 80 100 120 140 160 180 200 



Q [GeV] 

Figure 11: Tensor charge as a function of Q for up and down quark for the GRSV and GGR 
models. 
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